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Abstract 

We study the dependence of the heat transfer between two semi- 
infinite solids on the dielectric properties of the bodies. We show that 
the heat transfer at short separation between the solids may increase 
by many order of magnitude when the surfaces are covered by adsor- 
bates, or can support low-frequency surface plasmons. In this case the 
heat transfer is determined by resonant photon tunneling between ad- 
sorbate vibrational modes, or surface plasmon modes. We study the 
dependence of the heat flux between two metal surfaces on the elec- 
tron concentration using the non-local optic dielectric approach, and 
compare with the results obtained within local optic approximation. 

1 Introduction 

It is well known that for bodies separated by d » dx = ch/hsT, the 
radiative heat transfer between them is described by the Stefan- Bolzman 
law: 

S=^L(T?-T*), (1) 
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where 7\ and T 2 are the temperatures of solid 1 and 2, respectively. In this 
limiting case the heat transfer is connected with traveling electromagnetic 
waves radiated by the bodies, and does not depend on the separation d. For 
d < c/t the heat transfer increases by many order of magnitude due to the 
evanescent electromagnetic waves that decay exponentially into the vacuum; 
this is offen refered to as photon tunneling. At present there is an increasing 
number of investigations of heat transfer due to evanescent waves in con- 
nection with the scanning tunneling microscopy, and the scanning thermal 
microscopy (STM) under ultrahigh vacuum conditions [TJI21E1I11I3EIIZI- It is 
now possible to measure extremely small amounts of heat transfer into small 
volumes jSJ. STM can be used for local heating of the surface, resulting in 
local desorption or decomposition of molecular species, and this offer further 
possibilities for the STM to control local chemistry on a surface. 

The efficiency of the radiative heat transfer depends strongly on the di- 
electric properties of the media. In jSJ E] it was shown the heat flux can be 
greatly enhanced if the conductivities of the material is chosen to maximize 
the heat flow due to photon tunneling. At room temperature the heat flow is 
maximal at conductivities corresponding to semi- metals. In fact, only a thin 
film (~ lOA) of a high-resistivity material is needed to maximize the heat 
flux. Another enhancement mechanism of the radiative heat transfer can be 
connected with resonant photon tunneling between states localized on the 
different surface. Recently it was discovered that resonant photon tunneling 
between surface plasmon modes give rise to extraordinary enhancement of 
the optical transmission through sub- wavelength hole arrays |9J. The same 
surface modes enhancement can be expected for the radiative heat transfer 
(and the van der Waals friction JU]) if the frequency of these modes is suf- 
ficiently low to be excited by thermal radiation. At room temperature only 
the modes with frequencies below ~ 10 13 s _1 can be excited. For normal met- 
als surface plasmons have much too high frequencies; at thermal frequencies 
the dielectric function of normal metals becomes nearly purely imaginary, 
which exclude surface plasmon enhancement of the heat transfer for good 
conductors. However surface plasmons for semiconductors are characterized 
by much smaller frequencies and damping constants, and they can give an im- 
portant contribution to the heat transfer. Recently, enhancement of the heat 
transfer due to resonant photon tunneling between surface plasmon modes 
localized on the surfaces of the semiconductors was predicted in [7]. The 
authors studied the radiative heat transfer between a small particle, consid- 
ered as a point-like dipole, and a flat surface. However, this treatment can 



2 



be applied for scanning probe microscopy only in the case R >> d, where R 
and d are the radius of the particle and the separation between the particle 
and the surface, respectively. For the opposite limit R << d, which is more 
appropriate for the scanning probe microscopy, the heat transfer between the 
tip and the surface can, in the first approximation, be modeled by the heat 
transfer between two semi-infinite solids. In this case the multiple scatter- 
ing of the electromagnetic waves by the surfaces of the bodies, which was 
not taking into account in [7j, becomes important in the photon tunneling. 
In particular, at sufficiently small separation d, the photons goes back and 
forth several time in the vacuum gap, building up coherent constructive in- 
terference in the forward direction much as would occur in resonant electron 
tunneling. In this case the surface plasmons on the isolated surface combine 
to form a " surface plasmon molecule" , in much the same way as electronic 
states of isolated atoms combine to form molecular levels. This will result in 
a very weak distance dependence of the heat flux, because the transmission 
probability for photon does not depend on d in this case (see below). For 
large d the sequential tunneling is more likely to occur, where the photon 
excited in a surface plasmon mode, tunnels to the surface plasmon at the 
other surface, and then couples to the other excitations in the media and 
exit. Other surface modes which can be excited by thermal radiation are 
adsorbate vibrational modes. Especially for parallel vibrations these modes 
may have very low frequencies. 

All information about the long-range electromagnetic interaction between 
two non-contacting bodies is contained in the reflection factors of the elec- 
tromagnetic field. At present time very little is known about the reflection 
factors for large wave vectors and for extremely small frequencies. In our 
previous calculations of the radiative heat transfer and Van der Waals fric- 
tion [TTJ H21 El E] we mostly considered good conductors. In this case it was 
shown that the important contribution comes from the non-local optic effects 
in the surface region. However it was shown that the radiative heat transfer 
and Van der Waals friction becomes much larger for high resistivity material, 
for which the volume contribution from non-local effects is also important. 
Non-local optic refer to the fact that the current at point r depends on the 
electric field not only at point r, as it is assumed in the local optic approxi- 
mation, but also at points r' ^ r in a finite region around the point r. In the 
case when both points are located outside the surface region the dielectric 
response function can be expressed through the dielectric function appropri- 
ate for the semi-infinite electron gas. However, if one of the point r or r' is 
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located in the surface region, the dielectric response function will be different 
from its volume value, and this gives the surface contribution from nonlocal- 
ity. In order to verify the accuracy of the local optic approximation we study 
the dependence of the radiative heat transfer on the dielectric properties of 
the materials within the non-local dielectric approach, which was proposed 
some years ago for the investigation of the anomalous skin effects 



The problem of the radiative heat transfer between two flat surfaces was con- 
sidered some years ago by Polder and Van Hove pQ , Levin and Rytov j2] and 
more recently by Pendry |3] , and Volokitin and Persson E] ■ Polder and 
Van Hove were the first who obtained the correct formula for the heat trans- 
fer between two flat surface. In their investigation they used Rytov's theory 
[HI H3] of the fluctuating electromagnetic field. However, they presented 
their result only for identical media, and within the local optic approxima- 
tion. In the subsequent treatment, they made numerical calculations not of 
the heat flux itself, but its derivative with respect to temperature, i.e., their 
result pertain only to the small temperature differences. Unfortunately, their 
paper contains no analytical formulas in closed form . Levin and Rytov [2] 
used the generalized Kirchhoff's law [15 to obtain an expression for the 
radiative heat transfer between two semi-infinite media in the impedance ap- 
proximation. The case of the good conductors was investigated in the details 
both in normal and the anomalous skin effect region. Pendry [3| proposed a 
more compact derivation of the formula for the heat flux between two semi- 
infinite bodies due to evanescent waves, and calculated the heat transfer be- 
tween a point- dipole and a surface. Volokitin and Persson |2j considered the 
problem of the heat transfer between two flat surfaces, as a particular case of 
the general approach for calculation of the heat transfer. They investigated 
numerically the dependence of the heat flux on the dielectric properties of 
the bodies, and found that for good conductors, even for very small distances 
the heat flux is dominated by retardation effects. They also showed that the 
heat flux heat is maximal at conductivities, corresponding to the semi-metal. 

According to 121 EH 03 the heat transfer between two semi- infinite bodies 
, separated by a vacuum gap with the width d, is given by the formula 
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where 

1 fw/c 

+ 1 /7 x ^ + _ ,] (3) 

7H L/c | 1 — e 2 lPl d itip(u;)it2p(u;) p 

where the symbol [p — > s] stands for the terms which can are obtained from 
the first two terms by replacing the reflection factor R p for the p-polarized 
electromagnetic waves with the reflection factor R s for s- polarized electro- 
magnetic waves , and where p = ((cj/c) 2 — g 2 ) 1//2 , k — \p\. The Plank function 
of solid 1 

(w) = huj (e hul/kBTl -l)' 1 , (4) 

and similar for II2. The contributions to the heat transfer from the propa- 
gating (q < uj/c) and evanescent (q > uj/c) electromagnetic waves are deter- 
mined by the first and the second terms in Eq. (j3J), respectively. 

Let us firstly consider some general consequences of Eq. (jSJ). In the case 
oh heat transfer through free photons (q < uj/c), the transfer is maximal 
when both bodies are perfectly black and have zero reflection coefficient, 
R = R r + iRi = 0. Now, what is the photon- tunneling equivalent of a 
black body? For q > uj/c there are no constraints on the reflection coefficient 
R(q, uj) other than that lining, uj) is positive. Therefore, assuming identical 
surfaces, we are free to maximize the transmission coefficient corresponding 
to the photon tunneling 

D2 -2kd 

T = —thl (5) 

|1 -e~ 2kd R 2 \ 2 

This function is a maximum when [3] 

R 2 + R 2 = e 2kd (6) 

so that T = 1/4. Substituting these result in (j3J) gives the evanescent contri- 
bution 

k 2 T 2 n 2 

(C \evan _ rv B ± He /y\ 
\ D z)max ~ ^ ' 

where q c is a cut-off in q, determined by the properties of the material. It is 
clear that the largest possible q c ~ 1/a, where a is an interatomic distance. 
Thus, from Eq. © we get upper boundary for the radiative heat transfer at 
room temperature: (S z ) max ~ 10 12 Wm~ 2 . 
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We rewrite the denominator of the integrand in the term in Eq. (jSJ), 
which corresponding to the evanescent waves, in the form 

|1 - e~ 2kd R\ 2 = [(1 - e- kd R r ) 2 + e~ 2kd R 2 ] 

x [(1 + e- kd R r ) 2 + e~ 2kd R 2 } (8) 
The conditions for resonant photon tunneling are determined by equation 

Rr(u ± (q)) = ±e kd (9) 

This condition can be fulfilled even when exp(— 2kd) << 1 because for 
evanescent electromagnetic waves there is no restriction on the magnitude 
of real part or the modulus of R. This open up the possibility of resonant 
denominators for R 2 » 1. Close to resonance we can write 

(l±e- kd R r ) 2 ±e- 2kd R 2 

« e~ 2kd R' 2 {uj ± )[{uj -u ± ) 2 + {R l {u ± )/K^±))\ (10) 

where 



R r {u±) = 



UJ=UJ± 



which leads to the following contribution to the heat flux: 

The parameter q c in this expression defines the region < q < q c where 
the two- poles approximation is valid. To proceed further, let us make the 
following simplifications. Close to a pole we can use the approximation 

R = °— , (12) 

OJ — LUq — IT] 

where a is a constant. Then from the resonant condition Q we get 

lu± = lu ± ae~ kd . 

For the two poles approximation to be valid the difference Au = \uo + — u_\ 
must be greater than the width r\ of the resonance. From this condition we 
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get q c < hi(2a/r))/d. For short distances the parameter q c defines the value 
of q where the solution of Eq. (jUJ) ceases to exist. 
For lj > a and q c d > 1, from Eq. (fTTjl we get 

2 

J± = f £ [n 1 M-n 2 M]. (13) 

Interesting, the explicit d dependence has dropped out of Eq. (fT3|). How- 
ever, J may still be d- dependent, through the d- dependence of q c . For the 
small distances one can expect that q c is determined by the dielectric prop- 
erties of the material, and thus does not depend on d. In this case the heat 
transfer will be also distance independent. 



3 Numerical results 

Resonant photon tunneling enhancement of the heat transfer is possible for 
two semiconductor surfaces which can support low-frequency surface plasmon 
modes. The reflection factor R p for clean semiconductor surface is given by 
Fresnel's formula 



where 




where e is the bulk dielectric function. As an example, we consider two clean 
surfaces of silicon carbide (SiC). The optical properties of this material can 
be described using an oscillator model [TE] 

*M = (1 + / L ~/ T . v ) (16) 

with eoo = 6.7, to L = 1.8 • lO 14 ^ 1 , uj t = 1.49 • 10 14 s-\ and T = 8.9 • 10 U S - 1 . 
The frequency of surface plasmons is determined by condition e r (u p ) = — 1 
and from (JTSJl we get u> p = 1.78 • 10 14 s _1 . In Fig.l we plot the heat flux S(d): 
note that the heat flux between the two semiconductor surfaces is several 
order of magnitude larger than between two clean good conductor surfaces 
(see Fig. 3). 
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Another enhancement mechanism is connected with resonant photon tun- 
neling between adsorbate vibrational modes localized on different surfaces. 
As an example, let us consider ions with charge e* adsorbed on metal sur- 
faces. The reflection factor R p , which takes into account the contribution 
from an adsorbate layer, is given by [TTj : 

q-s/e + 47m a q[sa\\/e + qa±] 

tU, — : " ; : 7 , II/) 

q + s/e + Airn a q[sa\\/e — qa±\ 

and where au and a± are the polarizabilities of adsorbates in a direction par- 
allel and normal to the surface, respectively, e is the bulk dielectric function 
and n a is the concentration of adsorbates. For clean surfaces n a = 0, and in 
this case formula (J 17)) reduces to the Fresnel formula. The polarizability for 
ion vibration normal to the surface is given by 

e* 2 

a± = M(ui-^-iu Vl _y (18) 

where u±_ is the frequency of the normal adsorbate vibration, and rj± is the 
damping constant. In Eq. (fTTj) the contribution from parallel vibrations 
is reduced by the small factor 1/e. However, the contribution of parallel 
vibrations to the heat transfer can nevertheless be important due to the 
indirect interaction of the parallel adsorbate vibration with the electric field, 
via the metal conduction electron Thus, the small parallel component 
of the electric field will induce a strong electric current in the metal. The 
drag force between the electron flow and the adsorbates can induce adsorbate 
vibrations parallel to the surface. This gives the polarizability: 



a,, = — - ^ (19) 

n e (uo 2 — uj 2 — iur)\\) 



where n is the conduction electron concentration. As an illustration, in 
Fig. 2 we show the heat flux for the two Cu(001) surfaces covered by a low 
concentration of potassium atoms ( n a = 10 18 m~ 2 ) . In the q— integral in 
Eq.(jnj) we used the cut off q c ~ n/a (where a ~ Inm is the inter-adsorbate 
distance) because our microscopic approach is applicable only when the wave 
length of the electromagnetic field is larger than double average distance 
between the adsorbates. In comparison, the heat flux between two clean 
surface at separation d = lnm is two order of magnitude smaller. 
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Fig. 3 shows the thermal flux between two clean metal surfaces as a func- 
tion of electron density n. In the calculations we have assumed that one 
body is at zero temperature and the other at T = 273 K, and the Drude 
relaxation time r = 4 ■ 10~ 14 s _1 . When the electron density decreases there 
is transition from a degenerate electron gas (ksT << £p, where Ef is the 
Fermi energy) to a non-degenerate electron gas {k B T » e F ) at the density 
rip ~ (K B Tmf/ 2 /7T 2 h 3 , where m is the electron mass. At T = 273 K the 
transition density n F ~10 25 m -3 . The full line was obtained by interpolation 
between the two dashed curves, calculated in the non-local dielectric function 
formalism for the non-degenerate electron gas (valid for n < np ~ 10 25 m~ 3 ), 
and for the degenerate electron gas (for n > rip) [T^]. The thermal flux 
reaches the maximum S mSuX ~ 5 x 10 8 W-m~ 2 at n max m 10 25 m~ 3 , which cor- 
responds to the DC conductivity a ~ 3 • 10 3 (fi-m) _1 . Within the local optic 
approximation the radiative heat transfer is maximal at np maiX « 10 24 m -3 
where SL max ~ 10 9 W-m -2 . The thermal flux due to traveling electromag- 
netic waves is determined by formula which gives Sbb — 308 W-m -2 for 
T = 273 K. 

4 Summary 

We have studied the radiative heat transfer in dependence on the dielectric 
properties of the media. We have found that at sufficiently short distances 
between bodies the thermal flux can be significantly enhanced in compari- 
son with the black body radiation when the material involved support low- 
frequency adsorbate vibrational modes or surface plasmon modes, or the 
conductivity of the metals is chosen to optimize the heat transfer. This fact 
can be used in scanning probe microscopy for local heating and modification 
of surfaces. 
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FIGURE CAPTIONS 

Fig. 1. The heat flux between two semiconductor surfaces as a function 
of separation d . One body is at zero temperature and the other at T = 273 K 
For parameters corresponding to a surface of silicon carbide (SiC) (see text 
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for the explanation). (The log-function is with basis 10) 

Fig. 2. The heat flux between two surfaces covered by adsorbates , as 
a function of the separation d. One body is at zero temperature and the 
other at T = 273 K. For parameters corresponding to K/Cu(001) [TH] (u>± = 
1.9- lO 13 *- 1 ,^ = 4.5-10 12 s- 1 , 77|] = 2.8-10 10 s-\r/ ± = 1.6-10 12 s-\ e* = 0.88e) 
(The log-function is with basis 10) 

Fig. 3. The heat flux between two metal surfaces as a function of the 
free electron concentration n. One body is at zero temperature and the other 
at T = 273 K. The full line was obtained by interpolation between curves 
(dashed lines) calculated in the non-local dielectric formalism for a non- 
degenerate electron gas for n < rip ~ 10 25 m -3 , and for a degenerate electron 
gas for n > rip. Also shown is results (dashed lines) obtained within the local 
optic approximation. The calculation were performed with the damping 
constant r _1 = 2.5 • 10 13 s _1 , separation d = 10A and uq = 8.6 • 10 28 m~ 3 . 
(The log-function is with basis 10) 
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